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ABSTRACT 



Relying heavily on the Elsasser and Culbertson (1960) computa- 
tional system, computations of terrestrial radiation incident at a 
black body interface have been programmed for an arbitrary atmosphe- 
ric sounding. The program has been applied to a random sampling of 
62 radiosondes from the "Wark-sounding-catalog. " 

Since Wark et al. (1966) have kindly made available simulated 
Nimbus II channel 2 and channel 4 specific intensities for these same 
atmospheres, a multivariate regression was derived relating downward 
flux computations to Nimbus II channel 2 and channel 4 readouts as well 
as to two other gross air parameters: (a) total reduced water vapor 

depth, and (b) interface pressure. All four independent variables gave 
high statistical significance, with channel 4 filtered flux accounting 
for the major portion of the "explained variance" of the dependent 
variable. Total ozone was also tested but yielded no statistical sig- 
nificance on the regression. 
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1. Introduction 



A knowledge of downward infrared radiation flux is required in 
order to understand properly the heat partitioning in the earth-atmo- 
sphere and sun system. The atmosphere is much more transparent to 
short-wave than to long-wave radiation. This is depicted by the 
absorption spectra for the atmospheric gases on p. 153 of Fleagle and 
Businger [6] . From these spectra, one can see that there is strong 
absorption in the infrared by water vapor, carbon dioxide and to a 
lesser degree by ozone. It is the goal of this paper to develop a 
computer program that is readily adaptable to computation of downward 
radiation flux from soundings using input sounding data from each of 
these three atmospheric constituents. As a by-product, the secondary 
goal is to relate empirically the results to a regression equation of 
selected satellite air mass properties. 

Wark et al. [16] studied intensively the problem of upward flux 
and consequently compiled 106 randomly distributed radiosonde observa- 
tions from the I.G.Y. files in what is hereby termed the "Wark-catalog. " 
This composite set of soundings cover a wide range of radiative condi- 
tions. Moreover, simulated Nimbus II data was kindly provided by J.H. 
Lienesch of the National Satellite Center for the same 106 soundings. 
Downward infrared flux measurements may be made by the Kuhn-Suomi 
radiometer but these give only limited areal coverage compared to the 
MRR data from satellites; thus, it was decided to compute downward 
flux in a manner similar (not identical) to that employed by Wark 
et al. in his 1962 study, and to perform a multiple linear regression 
analysis with certain simulated Nimbus II channel -readout values as 
input variables . The computation of downward infrared flux is made 
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applicable to the top of the undercast cloud or to the ground level, 
but the computer programming is left flexible so as to apply to down- 
ward flux at an arbitrary number of sounding levels above the interface. 

2 . Nature of the Data 

The 106 radiosonde observations in the Wark-Gatalog are given as 
profile listings of pressure (mb) , temperature (°K) , mixing ratio of 
water vapor (g/kg ) * and mixing ratio of ozone (cm NTP per mb) . Each 
sounding or case atmosphere is identified by a number, the name of the 
radiosonde station, and the date and time of the sounding. The sky 
conditions is indicated by "clear" or by the height (in mb) of the 
top of the undercast cloud. Temperatures and pressures are given at 
approximately 25+5 selected levels as well as at the interface. 

The interface is defined for this study as being the earth's surface 
for clear cases or the top surface of the undercast for cloudy cases. 
Total (corrected) water vapor and ozone masses are also given in this 
Wark- catalog. The models developed by the Meterorological Satellite 
Laboratory were taken from radiosonde data over all latitudes and 
seasons. All the soundings were from the I.G.Y record-file and went 
to at least 25 mb. Above this level, temperatures were extrapolated 
parallel to the appropriate Supplementary Standard atmosphere to the 
pressure of 0.1 mb. Furthermore, stratospheric humidity was extrapo- 
lated to conform in general to the Supplementary Standard distribution 
for the particular air mass at the station (standard cP, standard mT 
air masses, etc.). The 62 cases selected are composed of 47 undercast 
cases at various altitudes and 15 clear cases, the latter being ran- 
domly selected from the Wark- catalog. 
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Nimbus II satellite data is used to compute filtered channel 4 
upward radiation flux and the interface Black Body radiation flux for 
each of the 62 selected atmospheres. The Nimbus II data lists total 
and filtered channel 4 (5y to 30y) infrared specific intensities for 
five zenith angles (theta), that is, 0=0°, 20°, 45°, 60°, and 78.5°. 

The specific intensities for 0 = 90° are obtained by a Lagrangian 
form of extrapolation which is explained below in connection with 
equations (19) and (20) . The radiation intensities of the 62 cases are 
listed in ergs/(cm 2 sec stdn) and are converted to watts/ (m 2 stdn) in 
this paper. It should be noted that the characteristics of the sat- 
ellite infrared sensing-filter system must be eliminated so that the 
results do not depend upon this optical system to any appreciable degree. 
The National Satellite Laboratory minimized the channel 2 and channel 
4 sensor-filter effects on the Nimbus II data by normalizing the filter- 
sensor response functions for the various I R wave channels, hence, 
the resulting response was the "best possible" estimate of the thermal 
flux (5-30y) for the channel 4 radiometer, and of the black body flux 
( 10 -lly) or its equivalent black body temperature for the channel 2 
radiometer. Also, the sounding data lists temperature at the inter- 
face of each of the case atmospheres, which Keith [9] found to be very 
accurately described by the Nimbus II channel 2 readout. 

3 . Computation of Downward Flux of Radiation 

Downward infrared flux of radiation arriving at the interface is 

quantitatively computed for each of the 62 selected atmospheres from 
the Wark-catalog by the Elsasser computational system. The downward 
flux computation is divided into five parts. 



15 



F d = + F C02 + F 03 - - AF (u* , \f ) (1) 

The last two terms of equation (1) represent the "so called overlap 
of flux" and, hence, are subtracted from the first three terms of 
equation ( 1 ) which represents the downward infrared flux not cor- 
rected in any way from the manner described by Elsasser and Culbert- 
son [5] (henceforth denoted by EC for abreviation) . The computation 
of the terms of equation ( 1 ) is explained below in connection with 
equations ( 6 ) through ( 11 ) . 

In actuality, an energy-sharing process occurs in the 15y and 
9.6y bands which correspond to the CO 2 -H 2 O and O 3 -H 2 O overlap fluxes, 
respectively, thus, accounting for the -AF-terms. For example, 
according to Hanel, Bandeen, and Conrath [7], ozone should be con- 
sidered the primary radiator in the 9.6y band, hence, water vapor 
is the weak secondary absorber-radiator resulting in the correction due 
to the shared radiation, Af (u* , y*) • 

It is necessary to know the layer values of reduced optical path 
for water vapor, ozone, and carbon dioxide in the computation of radia- 
tive fluxes by each of the three constituents. The reduced optical 
path du* of an absorber has been taken from EC, and its relation to the 
true optical path is given by 
du* = P ( ToV^du 

Po \ T / (2) 

According to EC [5], the dependence on temperature is slight and, hence, 
is excluded in this study. The finite difference relations for the 
resulting reduced optical paths of the three main constituents are given 
by 

K K 

u* = ZA-):u*i = E(q i+ }j) (Pi+Jj) (A i+ i«p) X 1.00705 X 10 - 6 (g/cm 2 ) (3) 

i=l i=l 
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(cm NTP) 



K 

U* = EA-jU* = 
i=l 



i=l 



K 

ECQi+ij) (A i+J5 p) X 10" 5 



(4) 



K 



K 



U* - EA/ll* = (^i+JjP) x 2.4469 X 10“4 (cm NTP) (5) 



Basically, the EC computational technique for the determination of 
downward infrared radiation flux is realized by the integration of 
the area under the curve as shown in Fig. 1 with the use of Elsasser- 
Culbertson [5] flux tables of R(u*,T) which have been incorporated 
in the author's computer program RADFLUX found in Appendix I-A. In 
essence, one enters the EC tables with logarithmic values of the 
constituent optical paths and their corresponding temperatures, and 
obtains the respective "R" values which are explained below in connec- 



tion with equations (6) through (11) . These "R" values are then 
integrated over the entire temperature range of the sounding result- 
ing in the flux of radiant energy due to the particular constituent. 
There are various forms that can be used to determine the energy 
integral but Elsasser et al. [5] and Yamamoto [17] used temperature 
instead of optical path (u) as the independent variable. 

The "R" term has been taken from EC and is given by 



hence, one can readily see from equation (6) that infrared flux can 
be realized by an integration scheme of the R(u*,T) term over the 
temperature range. For example, downward infrared radiation flux 
for water vapor and its corresponding finite difference form are 
given by 



i=l 



i=l 




( 6 ) 
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R(U,T) 



A 



Ref. 

Level 




< ,T (OK) 



Fig. 1. Schematic form of the distribution of R(U^T) from tables in 
Elsasser and Culbertson [ 5 ] as applied to a typical sounding representative 
of the U-T distribution for any of the three atmospheric constituents. 

The area bounded by the sounding and the isopleths U* = 0 and U* = Uf 
is the downward flux uncorrected for possible overlap. 
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(7) 



Tf T 

= S R(u f *,T)dT + ^ °r/u*(T),tJ dT 

o T f 

i=i f i=80°C 

F wv = 1 <V T i-l> + E Ri+l 5 ( u *f,T i )(T i+1 -T) + 

i=l i=i f 

f -273 o C 

-80C R(u f *,T)dT (8) 

For water vapor, the range of wave number (\b/ v f) is listed, by 40 
cm -1 intervals, from 20 to 2600 cm" 1 , which is regarded as spanning 
the effective infrared range of the terrestrial radiation spectrum. 
The finite difference approximation for the ^ in equation (8) 

assumes the form given by 



Ri + J* = [Ri(ui*,Ti) + Ri + i (ui + i*,Ti + x )] /2 



(9) 



and the last term of equation (8) is obtained from Table 20 of EC. A 
similar formulation for F C q 2 an<1 F 03 ls reallze< ^ by the use of Tables 
11 and 13, respectively, in EC; however, Rcc>2 and R03 val ues are re- 
stricted to the frequency ranges of 540-820 cm -1 and 970-1130 cm” 1 , 
respectively. 

The last two terms of equation (1) , the H2O-CO2 and H2O-O3 over- 
lap flux terms, are defined by the following expressions 

T f - T 0 _ 

AF(u*;U*) = f T F (u* f )R C02 (K* f ,T)dT + J t f (u*)R C0 (U*,T)dT (10) 

° T f 



AF(u*,U*)= J Tp(u* f ) R Q3 (U* f ,T)dT + J t f (u*) Rq 3 (U* ,T) dT (11) 



Moreover, values of Z^u*) in the H2O-CO2 overlap was determined by 
Martin [10] according to 



- v f v f ) 

Mu*) = jr Fv B vT Av = 



c= fL~^ 



-t2 



■'> 



dt) Av 



(12) 



^ B v ^Av 



vf 

^ B .t- Av 
v 0 vi 
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where L ^ is the generalized absorption coefficient listed by 40 cm“^ 
intevals and is the Planck specific black body flux for a given 

temperature . The results of this computation are shown in Table 1. 
The H 2 O-CO 2 overlap flux occurs in the 540-820 cm"^ range, hence the 
water vapor TF(u) within this range was obtained by performing the 
summation indicated by equation (12) using both T ^, and B ^ at 20°C, 
and at 40 cm"^ intervals. L v values were taken from Table 10 in EC. 
According to EC, 20°C is a representative base temperature for x 

Fv 

observations since the H 2 0 -C 02 overlap is in practice most important 
in the lower layers of the atmosphere. In contrast, Wark et al. [16] 
used Yamamoto's [17] Tp vo values which were based upon theoretical 
transition probabilities; however, Palmer* s [14] laboratory measure- 
ments made at 300°K seem to give good agreement with Yamamoto's 
values for 25 cm ” 1 intervals. After some additional smoothing, it 
was noted that the quoted L v values over 40 cm’^ intervals in Table 
10 of EC corresponded closely to Yamamoto's quoted values of "£ v /2" 
because of the slightly wider interval used in EC. 

Furthermore, ip(u*) in the H 2 O-O 3 overlap is given by 

Tp ( u* ) = exp [-0.1167 u*] (13) 

Equation (13) uses the continuum (970-1130 cm” 1 ) absorption coef- 
ficient of water vapor proposed by Hanel et al. [7]. In both the 
H 2 O-CO 2 and H 2 O-O 3 overlap regimes, the water-vapor sharing path has 
been increased by the factor 5/3 in order to transform beam trans- 
missivity to flux transmissivity. 

The results of the foregoing Elsasser computational scheme 
using the author's computer program RADFLUX on the 62 selected 
case atmospheres from the Wark-catalog are shown in Table 2 along 
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Table 1. Water vapor flux transmissivities in the H 2 O-CO 2 overlap 
regime (540-820 crn - ^-) . 



log U* Tf( u *) l°g u* xp(u*) 



-6.0 


0.9999997150 


-5.7 


0.9999994259 


-5.3 


0.999998527 


-5.0 


0.999997150 


-4.7 


0.999994259 


-4.3 


0.99998527 


-4.0 


0.99997150 


-3.7 


0.99994259 


-3.3 


0.9998527 


-3.0 


0.9997150 



-2.7 


0.9994259 


-2.3 


0.998527 


-2.0 


0.997150 


-1.7 


0.994259 


-1.3 


0.985665 


-.10 


0.97128 


-0.7 


0.94431 


-0.3 


0.86680 


0 


0.77328 


+0.3 


0.67567 


+0.7 


0.54833 
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with their corresponding case-number identifiers from this catalog. 
The table lists uncorrected fluxes for water vapor, ozone, and car- 
bon dioxide together with the CO 2 and ozone overlaps. The CO 2 and 
ozone overlaps are equally partitioned arbitrarily between CO 2 and 
water vapor, and ozone and water vapor, respectively. As a result, 
half of the overlapped flux is restored to CO 2 and ozone, that is, 
the partitioned overlap fluxes are subtracted from their respective 
uncorrected constituent fluxes as shown below in connection with 
equations (14) through (16), thus, realizing corrected constituent 
fluxes . 

Fwv(CORR) = F wv - AF(u*,U*)/2 - AF(u*,V)/2 (14) 

F C 02 (CORR) = F C02 - AF(u*,\C*)/2 (15) 

F 03 (CORR) = Fq 3 -AF(u*,U*)/2 (16) 

It should be noted that all the fluxes shown in Table 2 are in watts/ 
m and that the mean of each overlap and constituent flux as well 
as the average percent contribution of the corrected constituent 
flux on the total downward flux are given at the bottom of the table. 
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VALUE 


37 . A 2 3 


A. 1C7 


121.194 


luO.429 


8.779 


6.726 


43.780 


25.069 


132. 


223 


J»VG VALUe COM* IbUT ICN 


(Ji TOTAL 










2j. 07 






DOWNWARD FLUX 


IN PFRCEMT 


. 73.97 




5.96 









Table 2 . Selected 62 caa. atmospheres as lndantiflad from the Wark -catalog , and the computed 
downward flux from program RAD FLUX. 
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4 . Determination of Channel 4 Upward Flux of Radiation. 



For each sounding case, the filter-normalized upward flux in the 
wavelength of 5 y to 30y was available as a result of computations fur- 
nished by Lienesch. This was in the form of the "so-called" channel 
4 specific intensity at the top of the atmosphere. Channel 4 specific 
intensity I 4 ( 0 ) was then transformed to a channel 4 flux for each of 
the selected 62 cases by using the trapezoidal rule for integration on 

14 ( 0 ) as given by 

JL 2 

*F 4 = 2tt / 2 I 4 ( 0 ) d (sin 0 ) (17) 

o 2 

and its corresponding finite difference form is given by 



j TT / 2 

e i 4 (0 i )+i 4 (e i+ i) 

0i=O X%(sin 2 {0 i+1 } - sin 2 {0i>) 



(18) 



where i = 1, 2, 3,^, and 0 j. = 0°, 0 2 = 20 °, 0 3 = 45°, 0 4 = 60°, 



©5 = 78.5°, and 0^ = 90°. The trapezoidal rule was selected instead 
of other finite integration schemes because of the unequal zenith- 
angle intervals of integration. Also the channel 4 upward flux relation 
is the same as that used by Wark et al. which allows for the usual 

approximation of azimuthal symmetry. 

The zenith angles (0) of the Nimbus II listings went as far as 
78.5°, hence, the specific intensity f or 0 = 90° was extrapolated by a 
Lagrangian form expressed by 

1 ( 0 ) = Ii(0-0 2 ) (0-03) + I 2 (0-0]_) ( 0 - 03 ) + l3(0-0i) (0-0 2 ) (19) 

(01-02X01-03) ( 02 - 0 l) (02-03) ( 0 3 - 0 l) ( 03 - 63 ) 

and which reduces to 

I( 90 O) = I 78.50 + O. 6866 I 450 - 1.8649I 60 o + 1.1783l 785 o (20) 
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According to Wark et al. U6J. radiation for angles greater than 
78.5° contributes only about 4 percent of the total flux and the energy 
obtained by integrating eq. (17) between 0 = 78.5° and 0 = 90° lies 
well within this 4 percent range. 

The black body radiation flux (oTjJ) at the interface is given 
directly by the interface temperature (°K) , and when a= 0.56687 X 10“^ 
(watts/m^ °K^) , oT^ is in watts/m^. 

The two satellite air mass properties, <J>F^ and aT^, are made avail- 
able for this study by the foregoing procedures, where Tg is taken 
equivalent to the channel 2 emittance temperature in a simulated 
Nimbus readout. Moreover, the foregoing integration employing eqs. 

(18) and (20) were incorporated in the author’s computer program SUMFLUX 
(Appendix I C) . The 62 selected Nimbus II cases as well as the computed 
air mass properties are shown in Table 3. 

5 . Statistical Results and Inferences 

The object in the statistical treatment which follows is to obtain 
a specification equation for F^ in terms of satellite readouts and 
other readily accessible air-mass data. To accomplish this, there was 
available a BMD03R computer program that executes the multiple regres- 
sion analysis establishing a regression plane given by the linear form 
of 

Y = A q + A 1 Z 1 + A 2 Z 2 + A 3 Z 3 + A^Z^ 4* A 5 Z 5 (21) 

with the final results of the regression plane shown in Table 4. 

The original set of 62 independent data atmospheres was split 
randomly into two equal sets. One set, Sample A which is denoted by 
the asterisks in Table 3, was used to develop the regression while the 
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* * * * *• * * 



Table 3. Channel 4 filtered flux and interface black 
the 62 case simulated Nimbus II data. 

CHANNEL 4 
FLUX 



181.076 
186.376 
171.133 
153.014 
143.171 
189.664 
170.847 
163.286 
188.125 
155.662 
140.013 
105.440 
164.415 
156.219 
144.776 
163.762 
14C.085 
184.727 

120.077 
99.153 

133.449 

181.957 

141.883 

175.614 
162.747 
178.669 

96.914 

143.358 

145.276 

140.226 

106.511 

119.152 

122. 3 ei 

146.354 

151.586 

165.731 

176.863 
187.088 
164.452 
138.983 

131.614 
118.779 
155.021 

114.152 
101.118 
144.337 
143.031 
145.922 

116.863 
79.778 

117.937 
121.593 
130.516 
178.260 
175.544 
183.440 
167. 0C3 
180. 008 
174. 560 
163.219 
158.044 
132.669 



CASE 


INTERFACE 

PRESSURE 


INTERFACE 

TEMPERATURE 


* 2 


1CC0 


289 


* 3 


1009 


289 


* 4 


1014 


284 


* 7 


998 


267 


★ 8 


1000 


262 


* 10 


1000 


298 


* 12 


923 


282 


* 13 


1003 


278 


20 


850 


291 


23 


1000 


274 


27 


941 


258 


31 


1020 


235 


50 


850 


275 


51 


8 50 


268 


52 


500 


261 


53 


830 


275 


54 


7C3 


264 


55 


908 


285 


56 


526 


249 


57 


1 0C6 


226 


58 


400 


257 


59 


850 


291 


60 


500 


261 


61 


7C0 


283 


62 


652 


274 


63 


754 


283 


64 


400 


230 


65 


7C0 


267 


66 


700 


266 


67 


568 


261 


68 


500 


238 


69 


7C0 


248 



70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

92 

93 

94 

95 

96 

97 

98 

99 
*100 



700 
476 
466 
7C0 
810 
9 30 
932 
5 CO 
5C0 
4 CO 
8C6 
5C0 
400 
800 
722 
7C0 
7 20 
400 
818 
5C0 
3 70 
700 
7 CO 
850 
700 
850 
850 
724 
900 
700 



250 

264 

265 
280 
285 
289 
276 
256 
253 
244 
268 
242 
233 
264 
263 

266 
248 
221 
248 
246 

251 

284 
283 
289 
278 

285 
282 
273 
271 
255 



body flux on 



R LAC K 

BODY 

FLUX 

395.435 

395.435 

368.771 

288.090 

267.109 

447.042 

358.492 

338.581 

406.495 

319.511 

251.167 

172.884 

324.201 

292.430 
263.054 
324.201 
275.359 
373.993 
217.912 

147.883 
2^7.295 
406.495 
263.054 
363.604 
319.511 
363.604 
158.633 
288.090 
283.798 
263.054 

181.883 
214.432 

221.434 
275.359 
279.555 

348.430 
373.993 

395.435 
328.942 
243.469 
232.255 
200.929 

292.430 
194.422 
167.073 
275.359 
271.210 
283.798 
214.432 
135.224 
214.432 
207.598 
224.998 
368.771 
363.604 
395.435 
338.581 
373.993 
358.492 
314.872 
305.746 
239.687 
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remaining set, Sample B, was used to test the stability of the test 
data on the check-data which was also derived by the author’s computer 
program RADFLUX. 

According to Miller Cl3] and Martin et al. , the regression 

screening analysis uses a F’ critical value both to accept and remove 
individual variables from the regression. The critical F 1 value uti- 
lizes a ’’probability level” that is a function of the number ”k” of 
predictors tested and of the Sth selected prediction system. The pro- 
posed critical value corresponds to the tabulated F* value at a prob- 
ability level of a=a*/ (k-S-Fl) with 1 and N-S-l degrees of freedom, 
where N is the sample size to be tested. Since the predictors (k=5) 
have been pre-selected, the selection choice "S” is equal to 1, hence 
K-S+l is equal to 5 . 

The predictors employed are <f>F^. , oTg, u*, the interface pressure 
”p”, and U* which correspond to Zj_ 9 Z 2 , Z 3 , Z 4 and Z 5 , respectively, in 
eq. (21). The order of the predictors is a vital consideration in 
describing the ’’explained sum of squares by the regression.” Channel 
4 filtered flux (<^ 4 ) is referred to as the thermal channel, for its 
infrared range is 5y- 30y, whereas channel 2 filtered flux, the window 
channel, has an infrared range of 10-lly. Channel 2 filtered flux, 
however, is equivalent to the interface black body flux, oT^ ,and there- 
fore contains very little atmospheric structural information. Pre- 
scribing <j)F 4 ahead of oT^ allows <^4 to assume a larger proportion of 
the explained variance than oT^. 

The most readily adaptable test for the significance of the 
predictor is the F test based upon an analysis of variance. The F 
statistic is defined as 
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F = Mean squares explained by the predictor (s) (22) 

Mean squares not explained by the predictor (s) 

The F statistics derived from the results of the BMD03R program are 

shown in Tables 4 and 5. These statistics are used to determine the 

significant specification on F^ by the multivariate regression and by 

each of the predictors. Significance in this analysis means that the 

independent variable and/or variables selected for the F test did not 

affect the outcome of the regression by mere chance at the specified 

level of belief. For this analysis, an "a" of 0.01 is chosen resulting 

in a composite 99% confidence level of belief for the multiple regression. 

Also shown in Tables 4 and 5 are the tabulated critical F T values 
and the "a" critical level for rejection of a possible air-mass variable 
together with the corresponding "F-value upon entry." Since "a*" 
is chosen to be 0.01, the "a" critical levels for the 5 predictor 
regression are 0.002. The "F-values upon entry" are defined by 

F(l,N-m-l) = Mean Squares explained by the regression (23) 

Mean Squares not explained by the regression 

where "N" denotes the number of cases in the sample, "m" is the 

sequential order of the predictors (m = 1, 2....k), "k" is the number 

of predictors in the regression, and the number pair (l,N-m-l) denotes 

the number of degrees of freedom associated with the numerator and the 

denominator, respectively. For example, the "F-value upon entry" for 

the first predictor (^F^) in the regression is given by 

(24) 

F(l,29) = (Sum of Squares attributable to the regression)/! 

(Sum of Squares not attributable to the regression) /29 

The addition of the next predictor ,aTg , requires the definition of 

a new F-statistic specifically restricted to the testing of a Tg, 

after the contribution of ^F^ has been excluded. This new F-statistic 

follows the same form as given by eqs. (22) and (23) but with one less 
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degree of freedom in the denominator. The remaining independent 
variables are examined successively in the same manner with a new 
F-statistic specifically restricted to the testing of the selected 
variable after the contribution of the variables already tested have 
been excluded. 

The composite F-statistic for the multiple regression is based 

upon 

F(k,N-k-l) = (Sum of Squares attributable to the regression) /k [ ^ 

(Sum of Squares not attributable to the regression) / (N-k-1) 

Use of eq. (25) with the results from Table 4 gives the F-statistic 

for the 5 predictor multiple regression as F = 116.03 which is far in 

excess of a critical F ! (5,25) = 3.86 appropriate at a 99% confidence 

level . 

The results of Table 4 show not only that the multivariate 
regression, but also that each of the single variate predictors are 

significant with the exception of the last independent variable, U*. 

The "F-statistic upon entry" for u* is 0.90416 in contrast to the 
required (tabulated) critical F ! -value of 14.0 at 99.8% confidence level 
of belief. This indicates that U* could only have affected the outcome 
of the multivariate regression by chance at a very low level of probab- 
ility. An explanation for the lack of significance by ozone could lie 
in that the computed amounts of ozone in a column of atmosphere are near- 
ly constant across the 62 case sample resulting in a negligible influence 

on the regression. Since U* showed so little significance in the 
regression analysis, another analysis was made with U* excluded. The 
regression plane for the 4 predictor case is given by the linear 
expression 

Y' = Ai + A^Zx + A^Z 2 + A^Z 3 + A^Z 4 (26) 
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Table 4. Analysis of variance for the multiple regression in the five-predictor case. (Sample A) 
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116.031 3.86 0.9587 0.9791 -85.5347 



with the final results of the multivariate regression shown in 
Table 5. The same procedures and testing criteria as before were 
used resulting in each of the 4 predictors as well as the composite 
set of predictors showing significant contribution on the regression, 
thus, the 4 predictor regression was selected as the testing regres- 
sion plane for the independent data sample. The five variate regres- 
sion equation is given by 

F d = -85.535 + 1 . 750(4>F 4 ) - 0.486(aT^) + 47.466(u*) + 0.129(p) + 

-121.866(U*) (27) 

and the four variate regression equation is given by 

F d = -128.089 + 2.022(<t>F 4 ) - 0.605(oTj[) + 49.897(u*) + 0.129(p) (28) 
It is interesting to note that the multiple correlation coefficients 
shown in Tables 4 and 5 decrease very slightly from 0.9791 in the five 
variate case to 0.9784 in the four variate case. An explanation for 
this could lie according to Panofsky et al. [15] and Martin et al. in] 
that the longest equation may have actually overfitted the first analysis 
ascribing some of the variation due to small-scale perturbation to U* 
by chance. 

Finally, the 4 predictor regression was tested using Sample B 
data by an analysis of variance or of mean squares, that is, a test 
of variation about the regression plane. Again, significance of the 
regression plane was tested by the F test as given by equation (24) . 

Total sum of squares with its components is given by 
N _ N N 

z (y,-y) 2 = z (y.-y’.) 2 + z (y,-y) 2 (29) 

i=l i=i 1 1 i=l 

(A) (B) (C) 

where M y' " denotes the regression estimates for F d ,"yi M denotes the 
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Table 5. Analysis of variance for the multiple regression in the four predictor case. (Sample A) 



m 

CM 

o 

o 

o 

X II 

- 4-» 

B = 

a 

r— I 4-> “ 



G G 

u 



•H 


G 


> 








CM 




4-J 


G 


5-i 


no 


NO 




00 




•H 


i — 1 


4-J 


• 


• 


• 


• 




5-i 


G 


G 


O 


o 


o 


o 




O 


> 


G 


r — 1 


1— 1 


i—l 


1—1 






X 














G 


4-J 




CO 


o 


m 


CO 




G 


e 




00 


MT 


o 


o 




i — 1 




> 


m 


r-» 


00 


<r 




cd 


G 


U 


• 


• 


• 


• 




> 


O 


4-» 


o 


<1* 


CM 






i 


C4 


G 


o 


i — 1 


i—l 


CO 




Cn 


G 


G 


1—1 
















CM 


i — 1 




i—i 










ON 


O 


in 


CO 






C/3 




• 


• 


• 


• 




4-1 


G 


• 


NO 


i — 1 


00 


ON 




O 


5-i 


X) 


CM 


CO 


i—i 


NO 






G 


G 


O 


NO 




CM 




e 


G 


X> 




CM 


r^ 


ON 




G 


cr 


XI 


CM 


i—l 








CO 


GO 


c 


t — l 










1 — 1 






i— 1 


m 


m 


r- 




cd 


i— 1 


• 


ON 


00 








•H 


G 


4-1 


CM 


o 


m 


m 




4-» 


5-i 


4-1 




CO 








u 


5-1 


G 


• 


• 


• 


• 




G 


O 


O 


O 


o 


o 


o 






O 


O 




1 














NO 


1 — 1 


I—l 


CM 










« — 1 


m 




ON 










CM 


o 


ON 


CM 








• 


O 


NO 


00 


i—l 








44 


• 


• 


• 


• 




• 




4-1 


CM 


o 


ON 


O 




00 




G 






MT 






d) 




O 












Pd 




O 




i 














CO 


NO 


ON 


NO 


00 








ON 


o 


CM 


o 


m 








00 


NO 


ON 


NO 


00 


X> 




> 


r^. 


00 


o 


ON 


CO 


4-J 




G 


CM 


r-* 




o 




Crt 




Q 








CM 










NO 


nO 


NO 




in 








00 


NO 


00 




ON 








CO 




n-* 




O 




G 




• 


• 


• 


• 


• 




G 




ON 


CM 


o 


co 


CM 




G 




<d- 


ON 




<r 


co 




S 




i— i 


CM 






1—1 






G 












1 




i — 1 












5-t 




JO 












G 




G 


i — I 


CM 


CO 




— 


> 




•H 


CM 


CM 


CM 


CM 





4-> G 
G -H 
cd 

4-J g 

C/3 G 
cd i— i 
o cd 
o > 



c 

o 

•H 

C/3 

C/3 

G 

54 

00 

OJ 

Pd 



0) 

i — I • 

C4 ' — I 
*H G 
4J 54 
i— I 54 

G O 

S o 



<r 

Pd 



44 

44 

G 

O 

u 



C 

o 

•H 

03 

0) 

0) 

54 

00 

0) 

}-i 



0) 

X2 

4-J 



54 

o 

44 



03 

O 

•H 



03 

*H 

4-J 

cd 

4-J 

03 

4-J 

g 

0) 

G 

•H 

4-J 

u 

0) 

a. 

V-t 

0J 

X 

4-J 

o 



I 

4-1 cd Ov 
O G C\KT 
•H Pi 



44 5-1 

44 G G 
G 4-> O 
O 0) *H 

O Q 4-J 



i—l NO 

cd OJ CM 

•H i — I Ml* 
4-> g 
•H > 

V-i H - 

o P4 



C 



0) 
G 
« — I 
cd 
> 
I 

Cn 



no 

CM 

<r 



32 



145.3482 4.14 0.9572 0.9784 - 128.087 



observed values as determined by the computer program RAD FLUX, and "N ,r 
denotes the total number of cases in the independent data sample* Term 
(A) gives the total sum of squares; term (B) gives the sum of squares 
of deviations from the regression plane or the sum of squares not 
explained by the regression; and, term (C) gives the sum of squares of 
deviations of the regression estimates from the mean (y) of the observed 
values or the sum of squares explained by the regression. 

The statistics on the independent data sample are given in Table 
6 . The required (tabulated) critical F* value [>’ (4,26iJ is equal to 4.14 
at 99% confidence level of belief. One can readily see that the M F" 
statistic for the independent data F(4,26) = 118.2546, is considerably 
larger than its corresponding tabulated critical F f value, thus the 4 
predictor regression verifies as being significant for determining 
F^ within 99 % confidence level of belief. 

Furthermore, the results of the 4 predictor regression on the 
independent data sample were tested for shrinkage using the relation 
given by 

Rfc = 1 ~ (Sum of Squares not explained by the regression) (30) 

(Total Sum of Squares) 

* 2 

where is the effective coefficient of determination and R^ is termed 

the "effective multiple correlation coefficient." The results are shown 
in Table 6 . 

The percent shrinkage between the explained sum of squares of 

o - 2 2 

samples A and B is given by R 4 - R^ . Here is the coefficient of 

determination for eq. (28) (Sample A data). Table 6 shows that the 
regression performed exceptionally well on the independent data sample 
with only a 5.62% amount of variance shrinkage and with a R 4 and R 4 
of 90.1% and 95.1%, respectively. 
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Table 6 . Analysis of variance on the independent data (Sample B) 



Source of Variation 


Degrees 

of 

Freedom 


Sum 

of 

Squares 


Variance 


Explained Variance 


4 


86065.300 


2868.8430 


Unexplained Variance 


26 


9462.072 


157.6892 


Total Variance 


30 


95527.372 


3026.5322 





F-value 

F(4,26) 


Critical 
F* value 
F ! (4, 26) 


Coeff. of 
Determina- 
tion (R^) 


Multiple 
Correl. 
Coeff. (R 4 ) 


Percent 

Shrinkage 


118.2546 


4.14 


0.901 


0.951 


5.62 
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Table 7. Analysis of variance for the multiple regression in the four predictor case (Sample size 62) 
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270.0076 3.68 0.9499 0.9746 -129.847 



The foregoing analysis showed that the 4 predictor regression was 
indeed verified on the independent data sample; but, the process 
of randomly dividing the original 62 Nimbus II cases in order to realize 
this verification excluded 31 perfectly good cases from the regression 
analysis, thus, it was decided to obtain a regression based on the 
entire 62 Nimbus II cases. One should anticipate a more representative 
regression than the one previously derived, because the sample size is 
twice as large as that used in the original analysis. Although 
the new analysis used the same methods and testing criteria as that 
used in the 31 case sample, it was, however, confined to the 4 predictor 
analysis because of the earlier findings in the study. The final regres- 
sion is given by 

= -129.847 + 1.824(<t>F 4 ) - O^KaTg) + 43.293(u*) + 0.097(p) (31) 
and its corresponding statistics are given in Table 7. It is interest- 
ing to note that the multiple correlation coefficient decreased slightly 
in going from the 31 case analysis to the 62 case analysis. An explan- 
ation for this must lie in that the sample B "explained variance" was 
slightly smaller than that for sample A. 

6 . Conclusions 

The programmed computational system for downward infrared radiation 
flux due to the absorbing/emitting properties of water vapor, ozone and 
carbon dioxide is given in Appendix I A by the author's computer program 
RADFLUX. This program may be readily used for the determination of F^. 
The only necessary inputs are pressure (mb) , temperature (°K) , water 
vapor mixing ratio (g/kg) , and ozone mixing ratio (cm NTP per mb) for 
selected levels in a sounding. Moreover, this program can 1 be easily 
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modified as shown in Appendix I B in order to determine from any pre- 
selected level to the top of the atmosphere (0.1 mb). It should be 
noted that the mean ratio of downward radiation flux due to water vapor 

(F ) with and without overlap corrections to the total downward 
wv r 

radiation flux (F^) was 73.97% and 85.6%, respectively. 

Furthermore, total and filtered channel 4 upward radiation flux 
can be determined by using the author f s computer program SUMFLUX found 
in Appendix I C. Also, the interface black body flux or its channel 2 
equivalent temperature can be readily determined using this same pro- 
gram. 

Finally, the multivariate regression equation developed using 
Sample A Nimbus II air-mass properties for the determination of F^ 
had an accuracy of 95.1% on the independent test data. Sample B. This 
can be seen in the foregoing analyses. The explained standard devia- 
tion of the Sample A regression was 95.84%, but the resulting explained 
standard deviation on the independent sample was 95.1%, giving a loss 
of accuracy of 2.74%. A partial explanation for the loss of accuracy 
lies in the process itself of developing a regression on a finite ran- 
dom sample. One should not anticipate a regression determined from a 
finite sample to perform exactly the same way on another sample of the 
same population. This is inherent in the randomness of the selection 
process and of the variations within the population. This effect can 
be anticipated to be partly minimized by using the largest possible 
sample. Another partial explanation could lie in the assumption that 
the regression is linear, but in actuality, F^ might vary non-linear ly 
with one or any combination of its predictors. This would indicate 
that the true regression is some sort of surface and not a plane. 
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however, since the shrinkage of the multiple correlation was only 
2,74%, the regression plane appears to be a good estimate of the true 
surface or of physical reality. With the advent of improved satellite 
sensing systems, the regression equations can be expected to be 
improved by an analysis similar to the one given in this paper , 
Furthermore, an improvement on the regression plane may be anticipated, 
if the satellite data as well as the radiosonde data could be strati- 
fied to a common level, namely that of the interface. This, of course, 
would exclude pressure as a possible air-mass predictor, however, one 
could reasonably expect this stratification to reduce the scatter about 
the regression plane. Moreover, since the effect of ozone in a column 
of atmosphere gives a very nearly constant contribution to and, 
thus, has very little effect on the regression plane, and since one 
should expect that the regression based on the largest possible sample 
to be the most representative to the true regression surface, the 
best regression for the determination of F^ is, therefore, the 4 
predictor regression based upon the 62 case Nimbus sample, eq,»(31). 
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APPENDIX 1 



Significant computer programs developed for this study (Sections 
A, B, and C) . 
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Appendix I, Section A. Computer Program RAD FLUX (FORTRAN 63) 



COOP * BOX T. TUPAZ J . S/ 1 S/2S . 60 .40000 . 

FTN ,L »E« 

PROGRAM RADFLUX 

PURPOSE OF THIS STUDY IS TO EXAMINE DOWNWARD FLUX 
OD I MENS I ON TT< 13) .U(21 ) ,UU( 16 ) ,UUU( 23) .WATER! 2 1 . 13 ) .OZONE! 16.13). 
1C02 (23.13) .RWATER! 21 ) . ROZ ( 16 ) , RC02 ( 23 ) .P<40> . PMEAN ( 40 ) . PD I FF (40 ) . 

2 T ( 4 ) .Q!40) . QMEAN ( 40 ) . QDI FF ( 40 ) *CQ ( 40 ) .CQMEAN ( 40 ) .CQDIFF(40) . 

3SFW ( 40 ) .SFO! 40 ) .SFC(40 ) ,CFW<40) .CFO! 40) ,CFC( 40) .FUT! 40) .GUT (40) . 
4DUT ( 40 ) .FF ( 1 ) . GG ( 1 ) . DD ( 1 ) .COTAUF! 21 ) .03TAUF! 40) .RCOLAP! 1 ) . 

5C02LAP! 21 > .03RLAP(40) . R03LAP ( 1 ) ,UC (21) .COXMISI 21 ) .TCOLAP( 1 ) , 
6T03LAP ( 1 ) 

DIMENSION KASEI40) .LEVEL) 40) 

DIMENSION Z ( 40 ) . S ( 40 ) . V ( 40 ) ♦ R ( 40 ) 

OD I MENS I ON RH20< 40 ) .R03(40 ) ,RRCO(40 ) .TEMPI (40 ) ,TEMP2( 40 ) .TEMP3I40) , 
1R0PE! 13) . SOPE < 13) . TOPE ( 13) ,TTW( 13).TTO(13).TTC(13).RW(13).RO(13). 
2RC( 13) * RURF ( 13 ) .SURF ( 13 ) .TURF ( 13) 

DIMENSION MKASE ( 100) .ZCOLAPI 100) .Z03LAP! 100) .ZFWV! 100) .WFWV( 100) . 
1ZF03 (100) * WF03 ( 100) . ZFC02 ( 100) ,WFC02< 100) .ZTOTAL ( 100) 

DIMENSION RAT 101 ( 100 ) .RAT 1 02 1 100) .RAT 103 ( 100 ) 

READ! 50. 60 ) ( (WATER! I.J).J=1.13).I=1.21> 

READ! 50.60) ( ( C02 ( I.J) *J=1 *13) * 1 = 1*23 ) 

READ( 50.60) ( (OZONE! I . J ) . J= 1 . 13 ) . I = 1 . 16 ) 

60 FORMAT! 13F6.3) 

READ) 50.61 ) (U( I ) . 1=1 .21 ) 

READ! 50.61 ) (UUU! I).I=1.23) 

READ (50.61) (UU( I ) .1 = 1 .16) 

61 FORMAT! 12F4.0.3F3.0.F2.0.7F3.0) 

READ (50.59) (RWATER! K) ,K=1 .21 ) 

READ (50. 59) ( RC02 ( K ) . K.= 1 . 23 ) 

READ! 50.59) (ROZ(K) *K=1.16) 

59 FORMAT! 12F6.2 ) 

READ (50,701) ( COXM I S(I). 1=1.21) 

701 FORMAT ( 5F 15. 10) 

READ! 50.703 ) ( UC ( I ) , I = 1 . 21 ) 

703 FORMAT! 16F4.2 ) 

DO 62 1=1.13 
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62 TT ( I ) = I *10-90 
WRITE! 51,100) 

100 FORMAT! 1H1) 

WRITE! 51.63 )( TT( I). 1 = 1. 13) 

63 FORMAT (4X.13F8.0) 

WR I TE ( 5 1 . 100 ) 

WRITE! 51.64) ( U ( I ) .(WATER! I , J ) , J = 1 , 13 ) , I = 1 ,2 1 ) 

64 FORMAT (F5.1.13F8.3) 

WR I TE ( 51 . 100 ) 

WRITE! 51*64) < UU < I ) .( OZONE ( I . J ) , J=1 . 13 ) . I = 1 . 16 ) 

WRITE! 51.100) 

WRI TE! 51*64 ) (UUU! I ) * (C02! I , J ) , J=1 , 13 ) . I = 1 .23 ) 

INTER=63 
KUT = 0 
CONST = 0*4838 
777 READ (50.1000) KZ 
1000 FORMAT! 14) 

READ (50,7000) KOUNT 
7000 FORMAT! 14) 

10 READ (50,1) (KASE! I ) .LEVEL! I ) ,P( I ) ,T( I ) ,Q( I ) ,CQ( I ) . I =1 ,KZ ) 

1 FORMAT! 213, F9.1.F5.0.F12. 4, F10.3) 

WRITE ! 51 » 100 ) 

WRITE (51.2) (KASE! I ) .LEVEL! I ) ,P( I ) ,Q( I ) ,CQ( I ) * 1 = 1. KZ ) 

20F0RMAT! IX , 4HCASE , 4X . 5HLEVEL , 4X , 8HPRESSURE , 6X . 3HH20, 7X . 5H0Z0NE// ( IX 
1*14,18, F14.4.F12.4.F 11, 3) ) 

KT = KZ - 1 
SUMW=0. 

SUMO=0. 

SUMC = 0 • 

CONST1=0. 4342944819 
DO 20 1=1 »KT 

PMEAN! I ) = (P( I )+P( 1 + 1) )/2. 

QMEAN(I) =(Q( I )+Q( 1+1) )/2. 

CQMEAN ( I ) = (CO! I )+CQ( 1+1) )/2. 

PDIFF! I ) = P ( I ) -P ( 1 + 1) 

SlDIFF ( I ) = Q ( I ) -Q ( I + 1 ) 
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n r> n n 



CQD1F Ft I ) = CQ ( I )-CQ( 1+1) 

the purpose of thes sub-program is to calculate the log values for 
the reduced optical paths. 

5UMW=SUMW+QMEAN ( I )*PMEAN< I )*PDIFF! I )*( 1.00705E-6) 

SUMO=SUMO+CQMEAN! I )*PDIFF( I >*( l.E-5) 

SUMC=SUMC+PMEAN ( I )*PDIFF< I ) * ( 2 . 4469E-4 ) 

SFW< I ) =SUMW 
SFO ( I ) =SUMO 
SFC t I ) =SUMC 
IF (SFW(I)) 56.57.58 
58 CFW ( I ) =CONST l*LOGF ! SFW ( I ) ) 

56 CONTINUE 

IF ( SFO ( I ) ) 52.53.54 

54 CFO ( I ) “CONST l*LOGF ( SFO ( I ) ) 

52 CONTINUE 

IF (SFC (I)) 46.47.48 
48 CFC ( I ) “CONST l*LOGF l SFC (1)1 

46 CONTINUE 
GO TO 20 

57 CFW ( I ) *19. 

GO TO 56 

53 CFO ( I ) =19. 

GO TO 52 

47 CFC ( I ) = 19. 

20 CONTINUE 

OWRITEI 51.4) (I,T(I).P(I).0(I) , PMEAN ( I ) .PDIFF ( I ) .QMEAN ( I ) . S FW ( I ) . 

ISFOl I ) ♦ SFC ( I ) . I =1 »KT ) 

40FORMAT ( 2H I . 4X . IhT . 6X . 1HP . 7X . 1H0 »4X . 5HPMEAN » 3X . 5HPD I FF . 3X . 
15HC1MEAN.6X.3HSFW.6X. 3HSF0.6X.3HSFC// ( I 3 .F6 .0 . F8 . 1 »F 7 . 3 .2F8 . 1 *4F9. 4 
2 > ) 

WRITE! 51.100) 

WRITE! 51.9) ( I . T ( I + 1 ) .CFW! I ) »CFO( I ) .CFC! I ) . 1 = 1. XT ) 

90FORMAT 1 3H I .2X.6HT! I + 1 ) . 7X . 3HCFW . 17X . 3HCFO. 1 7X . 3HCFC/ / ! I5.F5.0. 
13E2 .6)) 

NOTE ANY LOG VALUE DENOTED AS 19 INDICATES THAT IT HAS BEEN 

FORCED OUT OF RANGE OF THE R VALUE TABLES BECAUSE OF THE U VA- 
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c lUE bEINo AERO. ALSO, NOTE THAT U VALUES CAN NOT BE NEGATIVE 

C DUE TO THE ARITHMETIC AND SEQUENTIAL PROCESSES OF THE PROGRAM. 

wRITEI 51,100) 

«JUMP= 1 
XH = 

NO = 

ICC = 

c the purpose of this sub-program is to determine THE SPECIFIC 

C WAT ER R VALUES FOR EACH LEVEL OF EVERY CASE. 

DO 13 I =KOUNT »KT 
TU = T( 1+1) - 273. 

*U=CFW! 1 ) 

3b DO 23 IC= 1 » 12 

IF(TU-TTdC) ) 23.40.50 
5 C I F ( TU-T T ( K+ 1 ) ) 40,43,23 
23 CONTINUE 

I F (JUMP )269.189»89 
S9 WR I TE ( 51 » 88 ) K.ASE ( I ) » T ( I + 1 ) »WU 

8 80 FORMAT ( IX » 8H CASE = I4.4X.8H TUH20 =F5.0,5X»11H LOG UH20 =F6.3.23 
1H ARE OUT OF TABLE RANGE) 

133 IF( JUMP! 107.26.13 
40 I T 

GO TO 65 
43 I T = K+ 1 
GO TO 65 

101 DO 39 L=1 .20 

IF(WU-UIL) ) 39.70.80 
80 IF ( WU-U ( L+ 1 ) ) 70.73,39 
39 CONTINUE 
GO TO 89 
70 IW = L 

74 I F ( JUMP >120,110,90 
73 IW = L+l 
GO TO 74 

900CALL FFINtWU.U! IW) ,U< IW + 1 > .TU.TT! IT > »TT ( IT + 1 ) .WATER I IW. IT ) . 

1 WATER! IW, I T+l ) .WATER! IW+1.1T), WATER! IW+1, IT + 1) »FUT< I ) ) 
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ICH = XHi-l 

RH2CHKH) = FUT ( I ) 

TEMPI ( KH) =T < 1 + 1 > 

13 CONTINUE 

WRITE! 51.100) 

JUMP=JUMP-1 
DO 26 I *ICOUNT *KT 
TU= T l I + 1 ) -273 • 

WU*=CFO( I 1 
GO TO 36 

189 WRITE! 51.188) (CASE ( 1 ) * T ( I + 1 > »WU 

1880F0RMAT! 1X.8H CASE = I4.4X.8F! TU03 = F5.0.5X.10H LOG U03 =F6.3,23H 
1ARE OUT OF TABLE RANGE) 

GO TO 133 

102 DO 105 L= 1 * 15 

I F t WU-UU < L ) >105*70.106 
106 I F ( WU-UU (L+l))70»73»105 
105 CONTINUE 
GO TO 189 

110LCALL FFOZtWU.UU! IW) ,UU( IW + 1 ) .TU.TT! IT > ,TTl IT + 1 ) .OZONE! IW. IT ) . 
1GZONE! IW. IT+1 ) .OZONE ( I W+l . IT ) .OZONE! IW+1 . IT+1 ) .GUT ( I ) ) 

XO=ICO+l 

R03 ! ICO ) =GUT ( I ) 

TEMP2 ( ICO > = T ! I + 1 > 

26 CONTINUE 
JUMP = JUMP- 1 
DO 107 I=KOUNT.KT 
TU= T ! I + 1 ) -273 • 

WU=CFCl I ) 

60 TO 36 

289 WR I TE ( 5 1 .288 ) RASE t I ) . T < I + 1 > .WU 

288 FORMAT! IX. 8H CASE = IA.4X.8H TUC02 = F5.0.5X.11H LOG UC02 =F6.3.23 
1H ARE OUT OF TABLE RANGE) 

GO TO 133 

103 DO 108 L= 1 . 22 

I F < WU-UUU ! L > > 108,70.109 
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109 IF (WU-UUU< L + l ) ) 70, 73, 108 
108 CONTINUE 
GO TO 289 

12 00CALL FFC02 ( WU »UUU ( 1 W > ,UUU< IW+1).TU*TT(IT),TT(IT+1),C02<IW.IT), 
1C02 I IW. IT + 1 ) .C02( IW+1, IT) .C02( IW + 1 , IT + 1 ) ,DUT ( I ) ) 

K.C=KC+1 

RRCOI <C ) = DUT ( I ) 

TEMP3IKC) = TII+l) 

107 CONTINUE 
GO TO 77 

6!> IF (JUMP )1G3»102»101 

this sub-program determines the flux integral for water vapor 

BETWEEN -80 TO -273 DEGREES CENTR I 6ADE • 

77 WU = CFWIKT) 

IND = 3 
DO 150 1=1,20 

IF (WU - U( I ) ) 150,170,180 

180 IF ( WU-U( 1+1 ) ) 170,173,150 

150 CONTINUE 

170 IW=I 

1 7 A GO TO ( 220,210, 190 ) IND 
173 IW=I+1 

GO TO 17A 

190 CALL RRWAT(WU,U( IW) ,U( IW+1 ) ,RWATER( IW) ,RWATER( IW+1) ,FF) 

WRITE (51,191) FF 

191 FORMAT ( IX ,6H FF = E17.9) 

THIS SUBPROGRAM DETERMINES THE FLUX INTEGRAL FOR OZONE BETWEEN 
-80 TO -273 DEGREES CENTRIGADE. 

IND=IND-1 
WU = CFO(ia ) 

BO 160 1=1,15 
I F ( WU-UU ( I) >160,170,181 

181 IF(WU-UU( I + 1 ) ) 170,173,160 
160 CONTINUE 

210 CALL RROZ(WU,UU( IW ) ,UU( IW+1 ) ,ROZ( IW) ,ROZ( IW+1 ) ,GG) 

WRITE (51,211) GG 
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211 FORMAT ( IX » 6H GG = E17.9) 

this sub-program determines the flux integral for co2 between 

-80 TO -273 DEGREES CENTR1GADE. 

1 ND = 1 ND-1 
WU = CFC< XT ) 

DO 200 1=1.22 
IF (WU-UUUI I ) ) 200,170,182 
182 IF (WU-UUUI 1+1 ) > 170,173,200 
200 CONTINUE 

220 CALL RRC02 ( WU *UUU ( 1 W ) , UUU ( IW + 1) ,RC02( IW) , RC02 ( I W+l ) , DD ) 

WRITE (51,221) DD 

221 FORMAT ( IX , 6H DD = E17.9) 

this su9-program integrates the integral rdt between t initial 

TO T FINAL IN DEGREES CENTRIGADE. 

1ND=IND-1 

areai=o. 

AREA2 = 0 • 

AREA3=0. 

total=o. 

XHH=XH-1 
DO 400 1*1, XHH 
Y*TEMP1( I ) - TEMPI ( 1+1 ) 

400 AREAl=AREAl+( RH20I I )+RH20( 1+1 ) ) /2.*Y 
WRITE (51,413) AREA1 

413 FORMAT ( IX, 10H AREA 1 = E17.9) 

XOO=XO-l 

DO 401 1*1, XOO 
Y = TEMP2 ( 1 ) - TEMP2 ( 1+1 ) 

401 AREA2=AREA2+(R03( I >+R03( 1+1) >/2.»Y 
WRITE (51,414) ARE A2 

414 FORMAT ( IX , 10H AREA2 = E17.9) 

XCC=XC-1 

DO 402 1*1, XCC 

Y*T EMP3 ( I ) - TEMP3 ( 1+1 ) 

402 AREA3=AREA3+( RRCOI 1 )+RRCO( 1 + 1 > ) / 2 • *Y 
WRITE (51,415) ARE A3 
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415 FORMAT) IX, 10H AREA3 = E17.9) 

c this sub-program integrates the integral rdt between t final 

C TO 80 OEGREES CENTRIGADE. 

MC = 

MO 3 
MW = 

KUMP*1 
WU=CFW ( KT ) 

00 900 L 3 1 » 1 3 
DO 903 J 3 1 » 20 
I F ( WU-U ( J ) ) 903,906,909 
909 IFIWU-U! J+l ))906,973,903 
903 CONTINUE 
GO TO 900 
906 IW = J 

974 IF(KUMP) 983,982,981 
973 I W= J+ 1 

GO TO 974 

981 CALL FFFN ( WU ,U ( IW) ,U( IW+1 > .WATER! I W, L ) » WATER < I W+l , L ) ,ROPE(L) ) 
MW=MW+ 1 

TTW ( MW > 3 TT ( L ) 

RW(MW) 3 ROPE ( L ) 

900 CONTINUE 

KUMP=KUMP-1 
WU=CFO( KT ) 

00 915 L= 1 , 1 3 
00 918 J 3 1 * 1 5 
I F ( WU-UU ( J ) ) 918,906,924 
924 I F ( WU-UU ( J+ 1 ) ) 906,973,918 
918 CONTINUE 
GO TO 915 

982 CALL FFFN ( WU.UU ( IW >.UU< I W + l) .OZONE! IW.L) .OZONE! I W + l ,L ) ,SOPE< L) ) 
M0=M0+1 

TTO(MO)=TT( L) 

RO(MO) = SOPE(L) 

915 CONTINUE 
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MJMP=KUMP-1 

*u=CFC(ia ) 

DO 930 L° 1 * 1 3 
DO 933 J«l. 22 
I F ( WU-UUU ( J ) ) 933 #906 t 936 
936 I F ( WU-UUU ( J+l ) ) 906*973*933 
933 CONTINUE 
GO TO 930 

98 3 CALL FFFN ( WU *UUU( I W ) *UUU( IW+1) *C02( IW.L) *C02 ( IW+1.L) *TOPE(L) ) 
MC=MC+1 

TTC(MC) = TT ( L ) 

RC(MC) = TOPE ( L ) 



930 CONTINUE 
MUMP=1 



TU = T EMP 1 ( KH ) -273. 

XL = ICH 
*M = XO 
♦CN = KC 

830 DO 800 £=1*12 

I F ( TU-TT ( K ) ) 800 • 803 • 606 
806 IF < TU-TT( £+1 ) ) 803*873*800 
800 CONTINUE 

I F ( MUMP ) 815*812*809 
809 WRITE! 51.818) TU 

818 FORMAT ( 10H TF H20 = F5.0.3X.22H IS OUT OF TABLE RANGE//) 
£L=£L~1 

TU=TEMPl(£L)-273. 

GO TO 830 
803 I T = £ 

874 I F ( MUMP ) 883*882*881 
873 I T = £ + 1 



GO TO 874 

881 CALL FORGE 1 (TU*TTW( IT ) *RW( IT) * RH20 ( £L ) .A RE All) 
ADD 1=0. 



I M = I T- 1 
DO 824 L= 1 * I M 
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CALL F0RGE2 (RW< L ) » RW ( L+l ) » T T W ( L ) * T TW ( L+ 1 ) »RURF ( L) ) 

82 A CONTINUE 

DO 827 L=1 * IM 
827 ADD1=ADD1+RURF(L) 

This sub-program determines the total flux for water vapor WITH- 
OUT C02 AND OZONE OVERLAP CORRECTIONS. 

845 AREA1 = AREA1 + ADD1 + AREA11 

WRITE (51.410) ADD1.AREA11 .AREAl 

410 FORMAT ( IX . 8H ADD1 - E17.9.5X.10H AREAl 1 = E17.9.5X. 14HT0TAL AREAl 
1= E17.9) 

MUMP=MUMP-1 
TU=TEMP2(KO) -273. 

GO TO 830 

812 WRITE) 51.833) TU 

833 FORMAT ( 10H TF 03 = F5.0.3X.22H IS OUT OF TABLE RANGE//) 

KM=KM-1 

TU=TEMP2< KM) -273. 

GO TO 830 

882 CALL FORGE l(TU.TTO< IT) .RO( IT) ,R03(KM) .AREA22 ) 

ADD2 = 0 . 

I M= I T- 1 
DO 836 L= 1 . IM 

CALL FORGE2 ( RO ( L ) .RO ( L+l ) . TTO< L) . T TO ( L+l ) .SURF ( L ) ) 

836 CONTINUE 

DO 839 L= 1 ♦ 1 M 
839 ADD2=ADD2+SURF( L) 

c This sub-program determines the total flux for ozone. 

848 AREA2'AREA2 + ADD2 +AREA22 

WRITE (51.411) ADD2.AREA22.AREA2 

411 F ORMAT ( IX »8H ADD2 = E17.9.5X.10H AREA22 = E17.9.5X.14H TOTAL AREA2 
1= E17.9) 

MUMP=MUMP- 1 
TU=TEMP3(KC) -273. 

GO TO 830 

815 WRITE( 51.842 ) TU 

842 FORMAT! 10H TF C02 = F5.0.3X.22H IS OUT OF TABLE RANGE) 
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AN=KN-1 

TU*TEMP3 ( K.N ) -273. 

GO TO 830 

883 CALL FORGE 1 (TU»TTC(IT)»RC(IT>. RRCO (K.N1.AREA33) 

ADD3=0. 

I M= I T-l 
DO 851 L=1 » IM 

CALL FORGE 2 ( RC ( L ) »RC(L+1 ) . TTC ( L ) »TTC ( L+ 1 » .TURF ( L ) ) 

851 CONTINUE 

DO 854 L=1 » I M 
854 ADD3=ADD3+TURF(L> 

c this sub-program determines the total flux for carbon dioxide. 

857 AREA3= AREA3+ADD3+AREA33 

rfRITE (51.412) ADD 3 * AREA33.AREA3 

412 FORMAT ( IX. 8H ADD3 = E17.9.5X.10H AREA33 = E17.9.5X.14H TOTAL ARE A3 
1= E17.9) 

C THIS SUB-PROGRAM DETERMINES THE C02 AND OZONE OVERLAP CORRECTIONS. 

JIP = 3 

DO 450 1=1.20 
WU=CFW ( KT ) 

I F ( WU-UC ( 1 I ) 450,470.480 
480 I F ( WU-UC ( I + 1 ) ) 470.473.450 

450 CONTINUE 

WRI TE< 51 .451 >WU 

451 FORMAT ( 17H LOG UF TAUC02 = F9.4.3X.22H IS OUT OF TABLE RANGE//) 
RCOLAP=0. 

GO TO 452 
470 I W= I 

474 GO TO (490.491.492 )JIP 

490 CONTINUE 

491 CONTINUE 
473 I W= I + 1 

60 TO 474 

4920CALL AOVLAP < WU ,UC ( IW ) .UCIIW+1 ) .COXMISl I W ) ,COXMIS( IW+1 ) .COTAUF(KT) . 
1 ADD 3 .AREA33.DD.RCOLAP) 

452 JIP=JIP-1 
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00 886 l=X0UNT.KT 

886 ©3TAUF ( I ) = EXPF(-0.1167 * SFwi I ) ) 

R03LAP=03TAUF (XT )* < ADD2+AREA22+GG) 

< I P = 3 
KV«= 

DO 549 J'XOUNT *KT 
WU = CFW< J) 

DO 550 1=1.20 
I F ( WU-UC ( I ) ) 550.570.580 
580 IF <WU-UC( 1 + 1 ) ) 570.573.550 
550 CONTINUE 

WRI TE( 51 .453 ) (CASE(J).WU 

4 5 30 FORMAT ( IX . 7H CASE =I4.3X.16H log UC02 TAU =» F9.4.22H IS OUT OF TAB 
IkE RANGE) 

60 TO 549 
570 IW = I 

574 GO TO (590.591.5921KIP 

590 CONTINUE 
573 IW=I+1 

GO TO 574 

5920CALL BOVLAP ( WU .UC ( IW) .UCl IW+1) ,COXMIS( IW) ,COXMIS( IW+1) . 

1C0TAUF( J) .OUT ( J ) .T ( J ) .T( J+l ) .C02LAPI J) ,DUT ( J+l ) ) 

KV=KV+1 

C02LAP ( XV ) = C02LAP ( J ) 

549 CONTINUE 

XIP = KIP - 1 
TCOLAP = RCOLAP 
DO 449 J= 1 * XV 

449 TCOLAP * TCOLAP + C02LAP(J) 

DO 454 J=X0UNT.XT 

591 CALL BBVLAP(03TAUF( J) .GUT( J) .GUT( J 1 ) » T ( J ) .T ( J+l ) .03RLAP ( J ) ) 

454 CONTINUE 

T03LAP * R03LAP 
DO 569 J*XOUNT .XT 
569 T03LAP = T03LAP + 03RLAP ( J ) 

c this sub-program determines the total downward flux for every 
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C SOUNDING. 

CALL ASUMF 1 ARE A1 .AREA2 .AREA3 .FF .GG .DD. TCOLAP .T03LAP . TOTAL ) 

WRI TE< 51 *79 ) ! AREA1 * AREA2 » AREA3 » FF *GG *DD * TCOLAP . T03LAP * TOTAL ) 
790FORMAT { 2X . 5HARE A1 . 3X . 5HAREA2 » 3X » 5HAREA3 . 6X .2 HFF . 7X . 2HGG. 6X . 2HDD . 
15X.6HTCOLAP*2X.6HTO3LAP.2X,10hTOTAL FLUX// ( 1X.7E 17.9 > ) 

FWV = AREA1 fFF-TC0LAP/2.-T03LAP/2. 

F03 = AREA2+GG-T03LAP/2. 

FC02 = AREA3+DD-TCOLAP/2. 

C THE FOLLOWING VALUES HAVE UNITS CAL/CM SQUARE/DAY 

WRITE (51.201) AREA1.FF*AREA2,GG*AREA3.DD*T COlAP * T03LAP .FWV.F03. 
1FC02 .TOTAL 

20 10FORMAT ( IX » 1 OH AREA 1 = E17.9//1X.1 H FF = E17.9//1X. 10H ARE A2 

1 = E17.9//1X.10H GG = E17.9//1X* 10H AREA3 * El 7 .9/ / IX . 10H DO 

2 = E17.9//1X.10H TCOLAP = E 1 7. 9/ / IX . 1 OH T03LAP = El 7. 9/ / IX . 1 OH 

3 FWV = E17.9//1X.10H F03 * El 7 . 9/ / IX . 10H FC02 = E17.9///1X. 

A23H TOTAL DOWNWARD FLUX = E17.9) 

XAREA1 = AREA1*C0NST 
XFF = FF*CONST 
XAREA2 = AREA2*C0NST 
XGG = GG*CONST 
XAREA3 = AREA3*CONST 
XDD = DD*CONST 
XCOLAP = T COLAP*CONST 
X03LAP = T03LAP*C0NST 
XFWV = FWV»CONST 
XF03 = F03*C0NST 
XFC02 = FC02*C0NST 
XTOTAL = TOTAL*CONST 

C The following values have units watts/meters square 

WR I TE ( 5 1 » 10 1 1 ) KASE(KT) 

1011 FORMAT ( IX. 7HCASE = IA) 

WRITE! 51* 10U1) XAREA1 .XFF* XAREA2.XGG*XAR EA3. XDD. XCOLAP. X03 LAP. XFWV 
1 .XF03.XFC02 .XTOTAL 

1001 FORMAT! IX. 10H NAREA1 = E17.9//1X.1 H NFF = E17.9//1X, 10H NAREA2 

1 = E17.9//1X.10H NGG = El 7. 9// IX . 10H NAREA3 = El 7.9// IX . 10H NDD 

2 = E17.9//1X.10H NCOLAP = E 1 7. 9// IX , 10H N03LAP = E17.9// IX. 10H 
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3NFWV = E17.9//1X. 10H NF03 » E 1 7. 9/ / IX . 10H NFC02 = E17.9//1X. 

423H TOTAL DOWNWARD FLUX = E17.9///I 
INTERS INTER-1 
KUT = K.UT + 1 
MKASE(KUT) = KASE(KT) 

ZCOLAP(KUT) = XCOLAP 
Z03LAP ( K.UT ) » X03LAP 

ZFWV(KUT) = XFWV 
ZF03 ( K.UT ) = XF03 
ZFC02IKUT) = XFC02 
ZTOTAL(KUT) = XTOTAL 
WFWVOCUT) = XAREA1 + XFF 
WF03 ( KUT ) = XAREA2 + XGG 
WFC02 ( KUT ) = XAREA3 + XDD 
IF (INTER) 888*999. 777 
999 WRITE (51.100) 

YKUT = KUT 
ADD1 = 0.0 
ADD2 = 0.0 
ADD3 = 0.0 



AVG11 


= 


0,0 








AVG22 


= 


0.0 








AVG33 


= 


0.0 








AVG44 


X 


0.0 








AVG55 


= 


0.0 








AVG66 


sz 


0.0 








AVG77 


= 


0.0 








AVG88 


= 


0.0 








AVG99 


= 


0.0 








DO 8050 


1=1 fKUT 






AVG11 


= 


AVG11 


+ 


ZCOLAPt I ) 


AVG22 


= 


AVG22 


+ 


Z03LAP ( I ) 


AVG33 


= 


AVG33 


+ 


ZFWV ( 


I ) 


AVG44 


= 


AVG44 


+ 


WFWV( 


I ) 


AVG55 


= 


AVG55 


+ 


ZF03 ( 


I ) 


AVG66 


= 


AVG66 


+ 


WF03 ( 


I ) 
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AVG77 * AVG77 + ZFC02U) 

AVG88 = AVG88 + WFC02(I) 

AVG99 = AVG99 + ZTOTAL(I) 

8050 CONTINUE 

AVG 1 = AVG11/YXUT 

AVG2 = AVG22/YXUT 

AVG3 = AVG33/YXUT 

AvG4 = AVG44/YKUT 

AVG5 = AVG55/YKUT 

AVG6 = AVG66/YKUT 

AVG7 = AVG77/YKUT 

AVG8 * AVG88/YXUT 

AVG9 = AVG99/YKUT 

00 1111 1*1. (CUT 

RAT 101(1) = ZFWV< I )/ZTOTAL( 1 > 

RA T I 02 ( I ( = ZF03I I )/ZTOTAL< I ) 

RAT 103 C I I * ZFC02( I ) /ZTOTAL ( I I 

1111 CONTINUE 

BO 1112 1 = 1 * ICUT 

AD01 = AD01 + RAT 101 ( I ) 

A002 = A0D2 + RATI 02 ( I ) 

ADD3 = ADD3 + RATI03I I ) 

1112 CONTINUE 

ADD 1 1 = t ADD1/YKUT >*100. 

A0D22 = (ADD2/YKUT 1*100. 

A0033 = ( ADD3/YKUT (*100. 

WRITE! 51.1022) ( MXASE ( I ) .ZCOLAP! I ) .Z03LAP! 1 ) .WFWV( I ) .ZFWV ( 1 ) . 

1 wF03 l I ) .ZF03 ( I ) .WFC02 t 1 ) .ZFC02( I ) .ZTOTAL! I I , 1=1 .(CUT ) 

1022 FORMAT ( 1X.85HCASE H20-C02 H20-03 UNCOR- CORRECT- UNCOR- CO 

1RRECT- UNCOR- CORRECT- T0TAL/1X.88H OVERLAP OVERLAP REC 

2 TED ED WV RECTED ED 03 RECTED ED CO 2 DOWNWARD/ 1 X . 8 

34h wv flux flux 03 flux flux C02 

4 Fi-UX FLUX FLUX / ( IX* I4.2F8.3.F10.3 .3F9. 3.3F10.3 ) ) 

WRITEI51.1031) AVG1 .AVG2.AVG3.AVG4.AVG5.AVG6.AVG7.AVG8.AVG9 
1031 FORMAT! 1 X . 4HMEAN/ 1 X . 5H VALUE . F7 . 3 . F8 . 3 . F 10 . 3 . 3 F9 . 3 . 3 F10 . 3 ) 
wRl TEt 51 . 1032 ) ADD11 .ADD22.ADD33 
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1032 FORMAT ( 1X*31HAVG VALUE CONTRIBUTION ON TOTAL/1X.32H DOWNWARD FLUX 

1 IN PERCENT .F6.2.F18.2.F20.2 ) 

688 END 

THIS SUBROUTINE EXECUTES FORWARD AND BACKWARD INTERPOLATION TO OBTAIN THE 
SPECIFIC ELASSER WATER R VALUE FOR EVERY REDUCED OPTICAL PATH/ TEMPERA- 
TURE PAIR. 

SUBROUTINE FFINIU.Ul .U2 . T .T 1 . T 2 »FU1 T 1 .FU1T2 ♦ FU2T 1 .FU2T2 .F ) 

FFUT1 = FU1T1+(FU1T2-FU1T1)*ABSF( I T-Tl ) / t T2-T1 ) ) 
FFUT2=FU2T1+(FU2T2-FU1T1)*ABSF< ( T-Tl ) / ( T2-T1 ) ) 

F1=FFUT1+(FFUT2-FFUT1) *ABSF ( (U-U 1 ) / (U2-U1 ) ) 

FFUT3 = FU1T2-(FU1T2-FU1T1)*ABSF((T2-T)/(T2-Tin 
FFUTA=FU2T2-IFU2T2-FU2T1)*ABSF< ( T2 T)/(T2-T1 ) ) 

F2= FFUT3+(FFUT4-FFUT3)*ABSF( ( U-Ul ) / ( U2-U1 ) ) 

F = ( FI + F21/2. 

END 

OSUBROUTINE FF0Z(U,U1 .U2 . T • T1 » T 2 » GU1T 1 .GUI T2 .GU2T 1 .GU2T2 »G) 

GGUT1 * GU1T1+(GU1T2-GU1T1)*ABSF((T-T1)/(T2-T1) ) 

GGUT2 = GU2T1+IGU2T2-GU2T1 )*ABSF( ( T-Tl )/(T2-Tl ) ) 
G1=GGUT1+{GGUT2-GGUT1)*ABSF( (U-Ul ) / (U2-U1 ) ) 
GGUT3=GU1T2-(GU1T2-GU1T1)»ABSF( (T2 T)/(T2-T1 ) ) 

GGUT4=GU2T2-( GU2T2-GU2T1 ) *ABSF ( (T2 T > / t T2-T1 ) ) 

G2 = GGUT3+ ( GGUT4-GGUT 3 ) *ABSF ( <U-U1 ) / (U2-U1 > ) 

G= ( G1+G2 ) /2 • 

END 

OSUBROUTINE FFC02 < U .U 1 ,U2 . T , T 1 , T2 » DU1T1 .DU1 T2 . DU2T 1 , DU2T2 . D > 

DDUT1 = DU1T1+(DU1T2-DU1T1 )*ABSF( ( T-Tl ) / ( T2-T 1 ) ) 

DDUT2 = DU2T1 + (DU2T2 - DU2T l ) * ABSF((T-T1) / (T2 - Tin 
D1 = DDUT 1+ ( DDUT 2-DDUT 1 ) *ABSF ( (U-Ul ) / (U2-U1 ) ) 

DDUT3 = DU1T 2-( DUlT 2 -DU IT 1 ) *ABSF ( (T2 T)/(T2-T1) ) 

DDUT 4= DU2T2-CDU2T2-DU2T1 )*ABSF( (T2-T )/ ( T2-T1 ) ) 

D2 =DDUT3 + ( DDUT 4-DDUT3 ) *ABSF ( ( U-Ul ) / ( U2-U1 ) ) 

D= ( D1+D2 ) /2 • 

END 

subroutine rroz ( u.ui .U2 .roui .rou2 »gg i 

GG = R0U1 + (R0U2-R0U1 )*ABSF< (U-Ul ) /(U2-U1 ) ) 

END 
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SUBROUTINE RRWAT ( U.U1 .U2 .RWU1 »R*U2 *FF ) 

FF=RWUl+< RWU2-RWU1 ) *ABSF< ( U-Ul ) / (U2-U1 > ) 

END 

Subroutine rrco2<u*ui »u2*rcui,rcu2 .ddi 

DD=RCU1+ ( RCU2-RCU1 ) *ABSF( ( U-Ul | / ( U2-U1 ) ) 

END 

subroutine asumf(a.b.c*d»e.f*x,y*z i 

2 = A+B+C+D+E+F-X-Y 
END 

SUBROUTINE BOVLAP ( U* U1 ,U2 * R1 »R2 *ZZ . D1 * T 1 . T2 * XX . D2 ) 
ZZ = R1 + (R2-R1)*A6SF( ( U-Ul ) / { U2-U1 ) I 
XX=ZZ*<Dl+02 )/2.*(Tl-T2) 

END 

SUBROUTINE AOVLAPIU.Ul »U2 . R 1 »R2 *Z . A . B * DD * X ) 

Z = R1 + ( R2-R 1 ) *ABSF ( (U-Ul > / ( U2— U 1 ) > 

X = Z * ( A+B+DD ) 

END 

SUBROUTINE BBVLAP < A* R1 *R2 *T1 *T2*Z ) 

Z= A* < ( R 1 + R2I/2. * (T1-T2H 
END 

SUBROUTINE FFFN < U »U 1 »U2 .R 1 *R2 , XX ) 
XX=R1+(R2-R1)*ABSF< ( U-Ul ) / ( U2-U1 I ) 

END 

SUBROUTINE FORGE1 ( TF . T 1 .R1 .RF , X > 
X=(R1+RF)/2.*(TF-T1> 

END 

SUBROUTINE FORGE2 ( R 1 . R2 » T 1 * T2 , Y I 
Y= ABSF ( (R1+R2)/2.*(T2-T1) » 

END 

END 

FINIS 

-execute. 



606 



58 



Appendix I B. Modification Instructions to Program RADFLUX (FORTRAN 63) 



1. Each sounding is subdivided into levels starting from the interface 
at level I = 1 to the top of the atmosphere (0.1 mb) at level I = KZ. 

KZ is the number of the last level in the sounding. 

I 

2. The following must be completed in order to determine F^ from any 
preselected reference level different from the interface level: 

a. Note the count or number of the pre-selected level in the sound- 
ing and introduce this number as shown below in the example. 

1. Suppose selected reference level is 7 out of 30 levels in 
a sounding, hence, punch this number on a data card with 
the unit digit in column 4. This reference level card is 
the second data card in each of the case atmospheres. Do 
this procedure anytime the reference level changes for any 
case-atmosphere . 

3. If a level is added or subtracted in a sounding listing, the number 
on the first data card of the respective case atmosphere must be changed 
accordingly. For example: 

a. Suppose that there are 29 levels in a sounding, thus, the first 
data card in the sounding must read 29; however, suppose that 1 
level is added to the sounding, thence, the first data card must 
be changed to read 30. 

b. One must place the unit digit of the specific number in column 
4 of the data card. 

4. Finally, the INTER (constant) card must replaced when the number 
of case atmospheres in the study changes. For example: 

a. Presently there are 63 case atmospheres as data input to program 
RADFLUX, thus, the INTER card reads "INTER = 63." 
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b. Suppose the data input changes to 100 case atmospheres; thence, 
the INTER card must be changed to read "INTER = 100." 

5. Program RADFLUX was written for FORTRAN 63, however, this program 
can be made operational in any FORTRAN computer language provided 
that minor modifications to the control cards are made. One must 
check with the respective computer language manuals for further 
information. 
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Appendix I, Section C. Computer Program SUHFLUX (FORTRAN 63) 



-COOP , BOX T » TUPAZ J *S/ 1S/2S .20 . 20000 . 

-FTN.L.E. 

PROGRAM SUMFLUX 

The purpose of this program Is to integrate radiation intensity 

CASE. 

BETWEEN 0 AND 90 DEGREES TO OBTAIN TOTAL UPWARD FLUX FOR EVERY 
ODI MENS I ON KASE! 7) .PRES (7 ) .THETA (7) .X(7) ,X2(7 ) ,UH2 (7) »U03 ! 7 ) ,X4 ( 7 ) , 
1TEMPI7) .TOT INTI 7) ,CH4INT(7) 

DIMENSION KOUNT! 100) .PRESS! 100) .TEMPER ( 100 ) .CHFLUX! 100) . 

1BBFLUXI 100) 

INTER = 62 
HUT = 0 

CONST=0 • 56687E-7 
PI * 3.14159 
RADCON = 0.017453 
KZ = 6 

THETA( 6 ) = 90.0 

7770READ ( SO .2 ) KASE< 1 ) . PRES ( 1 ) . THETA ( 1 ) »X ( 1 ) »X2 ( 1 ) .UH2t 1 ) .U03( 1 ) »X4( 1 
1 I .TEMPI 1) 

2 FORMAT! 1 3 .F9 . 1 . F8 . 0 ♦ F 10 .3 .2F7 . 3 .F6 . 3 .F 10 . 3 .F7 . 1 ) 

K.T = KZ-1 

OREAD (50,3) < KASE <I) » THET A < I ) , X (I) »X2 ! I ) ,X4 ( I ) , I = 2 .KT ) 

3 FORMAT ( I3»9X»F8.0»F10.3»F7.3»13X»F10.3 ) 

SUM 1=0. 

SUM2=0 • 

DO 100 1=1, KT 

CALL FORD1 1 XI 3) »X ( A ) . X ( 5 ) »X4 ( 3 ) »X4 ( 4 ) .X4 ( 5 ) »X ( 6 ) , X4 ( 6 ) ) 

CALL FORD2I XI I)»X(I+1)»X4(I>»X4(I+1)»THETA(I),THETA(I+1).T0TINT<I) 
1 ♦CHAINTU ) .RADCON) 

WRITE! 51 .5 ) KASE ( I ) ♦ I .TOT I NT ( I ) .CHAINT ( I ) 

5UFORMAT! 1X.8H CASE = I5.5X.5H I = IA.5X.10H TOTINT = E17.9.5X.10H C 
1 HA I NT = E17.9//) 

SUM1=SUM1+T0T INT ! I ) 

SUM2=SUM2+CH4INT! I ) 

100 CONTINUE 

WRITE! 5 1.4) (KASE! 1) .PRES! 1 ) .TEMP! 1 ) . < I .THETA! I ) » X < I ) ,XA < I ) , I = 1 ,KZ 
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1 ) ) 

40F0RMAT < 1X.8H CASE = I5.5X.12H PRESSURE * F5.0.5X.8H TEMP = F5.0// 
1<1X,5H I = I4.5X.10H THETA1 * F5.1.5X.19H TOTAL INTENSITY = F6.3. 
25X.18H CH A INTENSITY = F6.3//)) 

FLUX1«SUM1*PI 
F LUX 2 “SUM 2 *P I 

WR I TE < 5 1 .6 ) RASE! 1 > .FLUX1 .FLUX2 

6 FORMAT! IX. 8H CASE * I5.5X.9H FLUX1 = E17.9.5X.9H FLUX2 = E17.9///) 
CALL FORD 3 ( CONST .TEMPI 1 ) .FLUX 3 ) 

WRITEI51.300) FLUX 3 
300 FORMAT! IX. 9H FLUX3 = E17.9///) 

PRESLN = LOGFIPRES! 1) ) 

INTER = INTER-1 
RUT = RUT + 1 
ROUNT < RUT ) = RASE! 1) 

PRESS! RUT ) = PRES! 1 ) 

TEMPER (RUT) « TEMP(l) 

CHFLUX! RUT ) = FLUX2 
BBFLUX! RUT ) = FLUX3 
IF (INTER) 888.999.777 
999 WRITE! 51.200) 

200 FORMAT ( 1H1 ) 

WRITE! 51 .400 ) ( ROUNT ( I ) .PRESS! I ) .TEMPER! I ) .CHFLUX! I ) .BBFLUX! I ) . 1 = 1 

1 .RUT ) 

400 FORMAT! IX. 53HCASE INTERFACE INTERFACE CHANNEL 4 BLACR 

1/1X.52H PRESSURE TEMPERATURE FLUX BODY /49X » 

2AHFLUX//I 1X.I4.F8.0.F13.0.F15.3.F14.3) ) 

888 END 

SUBROUTINE F0RD1 ( A .8 »C » D. E .F »G . H ) 

G = C+0.6866*A -1.8649*B +1.1783*C 
H = F+0.6866*D -1.8649*E +1.1783*F 
END 

SUBROUTINE FORD 2 (A.B.C.D.E.F.X.Y.Z ) 

X=(A+B) / 2 * * ( <SINF(F*Z ) )**2 - ( S I NF < E*Z ) ) **2 ) 

Y = ( C+D) /2.* ( ( SINF! F*Z ) ) **2 - < S I NF ( E*Z » ) ** 2 ) 

END 
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SUBROUTINE FORD3(A.B,C> 

C»A«B*»4 

END 

END 

FINIS 



-execute 
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